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Abstract. The amplitude for a spin foam in the Barrett— Crane model of Riemann- 
ian quantum gravity is given as a product over its vertices, edges and faces, with one 
factor of the Riemannian lOj symbols appearing for each vertex, and simpler factors 
for the edges and faces. We prove that these amplitudes are always nonnegative for 
closed spin foams. As a corollary, all open spin foams going between a fixed pair 
of spin networks have real amplitudes of the same sign. This means one can use 
the Metropolis algorithm to compute expectation values of observables in the Rie- 
mannian Barrett-Crane model, as in statistical mechanics, even though this theory 
is based on a real-time (e'"^) rather than imaginary-time (e""^) path integral. Our 
proof uses the fact that when the Riemannian lOj symbols are nonzero, their sign 
is positive or negative depending on whether the sum of the ten spins is an integer 
or half-integer. For the product of lOj symbols appearing in the amplitude for a 
closed spin foam, these signs cancel. We conclude with some numerical evidence 
suggesting that the Lorentzian lOj symbols are always nonnegative, which would 
imply similar results for the Lorentzian Barrett— Crane model. 



Physicists have long found it difficult to deal with the oscillatory behavior of the 
integrand in real-time Feynman path integrals of the form 



For this reason, we have become accustomed to using Wick rotation to switch to imag- 
inary time, thus replacing these integrals by those of the form 



where the positivity of the exponential allows one to use concepts from probability 
theory. This is helpful not only in analytical work but also in numerical computations. 
For example, computations in lattice field theory rely heavily on techniques such as the 
Metropolis algorithm. They give powerful tools for computing the latter sort of integral, 
but are typically useless for those of the former sort. 

Unfortunately, the strategy of Wick rotation becomes problematic in the context of 
quantum gravity. The basic difficulty is that in the absence of a favored time coordinate, 
the replacement t n- i< makes no sense unless one can prove that the results obtained 
thereby are independent of the choice of the t coordinate. A great deal of work on 
'Euclidean quantum gravity' is of dubious relevance to real-world physics because of 
inattention to this issue. 

Luckily, in the absence of topology change (i.e. for spacetimes of the topology K. x S") 
there exist theorems justifying Wick rotation for certain diffeomorphism-invariant field 
theories . Taking the prohibition against topology change seriously has led to models 
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of quantum gravity in which Wick rotation is indeed justified, and these appear to be 
better-behaved than those previously considered in EucUdean quantum gravity |^ . 

Nonetheless, there are still reasons to be interested in real-time Feynman path inte- 
grals in quantum gravity. After all, the prohibition against topology change might be a 
limitation of the Wick rotation method rather than an actual law of physics. Moreover, 
the original Feynman integrals have a certain conceptual priority over their Wick-rotated 
counterparts: quantum theory is ultimately about amplitudes, not probabilities. 

Yet another reason is that there are some very interesting models of quantum gravity 
based on real-time path integrals: namely, spin foam models ||, |lO|, |2^, p3[ . 

In a spin foam model of quantum gravity, the real-time path integral is rewritten as a 
sum over 'spin foams': 2-dimensional cell complexes with faces labelled by group repre- 
sentations and edges labelled by intertwiners. This sum has been rigorously shown to 
converge in some of these models at least for a fixed triangulation of spacetime. 

The next step is to determine whether these models have good behavior at large distance 
scales: that is, whether they reduce to general relativity in a suitable limit. For this, 
we need to compute expectation values of observables, e.g. by numerical methods. 

Here we show that for a specific model of this sort — the Riemannian Barrett-Crane 
model — one can do such computations using the Metropolis algorithm, because the 
amplitudes for closed spin foams are nonnegative. This model is not physically realis- 
tic: it is a version of 'Riemannian quantum gravity', that is, the quantum field theory 
based on a real-time path integral involving the Einstein-Hilbert action for Riemannian 
metrics. Since the Cauchy problem is ill-posed in classical Riemannian general rela- 
tivity, we should expect to run into problems with the quantum theory at some point. 
Nonetheless, the Riemannian Barrett-Crane model has so far proved to be a useful 
guide to work on the physically more important but mathematically more challenging 
Lorentzian Barrett-Crane model. In particular, we conjecture that the nonnegativity 
of spin foam amplitudes holds for the Lorentzian model as well as the Riemannian one. 

In Section 2 we prove the nonnegativity of spin foam amplitudes in two versions of 
the Riemannian Barrett-Crane model: the version due to De Pietri, Freidel, Krasnov 
and Rovelli pq| , and the version with modified edge and face amplitudes due to Perez 
and Rovelli [22[. In the Perez-Rovelli version the partition function has been shown 
to converge |2l]; in the DFKR version the partition function appears to diverge 
However, in both versions the spin foam vertex amplitudes are given by the Riemannian 
lOj symbols, and the key fact we need to show is that the product of these lOj symbols 
over all vertices of a closed spin foam is nonnegative. We prove this in two steps. In 
Lemma ^, we show that when the Riemannian 10 j symbol is nonzero, its sign is positive 
or negative depending on whether the sum of the ten spins in question is an integer or 
half-integer. In Lemma |^, we show that all the resulting minus signs cancel when we 
take the product of lOj symbols over all vertices of a closed spin foam. This immediately 
gives the desired result. As a corollary of this result, we show that all open spin foams 
going between fixed spin networks have real amplitudes of the same sign. Finally, we 
present some numerical evidence supporting our conjecture that the Lorentzian lOj 
symbols are always nonnegative. 
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2. The Riemannian Model 
In what follows we shall freely make use of standard graphical notation for SU(2) 



spin networks [13, 17 . We normalize trivalent vertices so that 



whenever the spins ji,j2,i3 form an admissible triple: that is, satisfy the triangle in- 
equality and sum to an integer. (When the spins do not form an admissible triple, we 
define the trivalent vertex to be zero.) As this normalization sometimes requires divid- 
ing by the square root of a negative number, it introduces a potential sign ambiguity. 
Luckily, in our calculations trivalent vertices always come in matching pairs, so these 
signs cancel, leaving an unambiguous result. This is important, because signs are of the 
essence in all to come. 

We shall also use the standard notation for simple spin networks ||l^, In such 
a spin network, labelling an edge by the spin j really means that it is labelled by irre- 
ducible representation j(8)j of the group Spin(4) = SU(2) x SU(2). Such representations 
are called 'simple', and they all descend to representations of S0(4). Unlike SU(2) spin 
networks, which change sign when one inserts a twist in the framing of any edge la- 
belled by a half-integer spin, simple spin networks are framing- independent. In a simple 
spin network, any unlabelled 4-valent vertex is really labelled by the Barrett-Crane 
intertwiner. This intertwiner is defined in terms of SU(2) intertwiners by: 





E 



(-l)2^(2j + l) 





where we sum over spins j such that the triples ji,j2,j and js, j4, j are both admissible. 
Superficially, this definition involves splitting the four edges incident to the vertex into 
a pair of 'incoming' edges and a pair of 'outgoing' ones. However, the Barrett-Crane 
intertwiner is actually independent of this choice. For example, we also have 




E 



(-lf^-(2j + l) 




a formal proof was given 



This fact is stated in Barrett and Crane's original paper 
by Yetter |6|. 

It will be technically useful to modify the Barrett-Crane intertwiner by removing the 
signs from the above formulas. Our final results will concern the original Barrett-Crane 
intertwiner, but our intermediate calculations will use the intertwiner without signs. In 
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this modified intertwiner we use a dotted line to indicate the sphtting of incident edges 
into 'incoming' and 'outgoing' pairs: 



In fact, this modified intertwiner only differs from the usual one by an overall sign. 
Lemma 1. 



(_l)2(ii+j2) 



(^_l)2(j3+j4) 



It follows that 



unless ji + ■ ■ ■ + ji is an integer. 

Proof — Since the triples ji,j2,j and are admissible in the sum 



,31 32 



31 32 



31 32 



^(-ir^(2j + l) 



we have (-1)^-' = = (_i)20'3+j4)^ which yields the equations to be proved. 

□ 

This means that the modified intertwiner changes at most by a sign when we change 
the splitting. For example, consider 



,31 32 



33 34 




= ^(2j + 1) ).^.: 









34 



This satisfies: 
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Lemma 2. 




J3 J4 



J3 Ji 



33 J4 



Proof — By Lemma |^ we have 



(-1) 



2(ii+j2) 



J3 J4 



J3 J4 



Rotating the pictures clockwise and relabelhng, we have 



^2(i2+J4) 



33 Ji 



33 34 



Combining these equations, it follows that 



Jl 32 / \3l 32 ^ 



33 34 



33 34 



as desired. A similar argument where we rotate counterclockwise gives 



31 32 



31 32 



(-1) 



2(J2+J3) 





□ 



The usual Riemannian lOj symbol [g^ is the function of 10 spins obtained by evalu- 
ating the following closed SU(2) x SU(2) spin network: 
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where the vertices are labehed by Barrett-Crane intertwiners. Mimicking this con- 
struction, we can also define a modified lOj symbol using the modified Barrett-Crane 
intertwiners: 




Here the choice of splittings really matters: by Lemma ^, if we change the splitting at 
any vertex, the value of the spin network might change sign. The nice thing about the 
above choice is that it always gives a nonnegative result: 

Lemma 3. 




> 



Proof — Using the definition of the modified Barrett-Crane intertwiner, the modified 
lOj symbol above can be rewritten as the following weighted sum of SU(2) spin networks: 

k2 k2 



5 k 

j2 (n(2fc.+i)) 



U J4 

Since the value of a disjoint union of spin networks is the product of their values, it 
follows that 




(1) 




= E (17(2^^ + 1)) 
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To show that the modified lOj symbols are positive, it thus sufHces to show that this 
SU(2) spin network, caUed a 15j symboh 



evaluates to a real number. 

To see this, first note that the only place where complex numbers are essential to the 
theory of trivalent SU(2) spin networks is in the normalization of the trivalent vertex. If 
we used the projection from the tensor product ji ®i2 to the summand j'3 as our trivalent 
vertex whenever ji, j2, js is an admissible triple, all closed trivalent SU(2) spin networks 
would evaluate to real numbers, since we could define them using the representation 
theory of SL(2,]R) instead of SU(2). However, this would not give 



Writing 0{ji, j2, js) for the unnormalized evaluation of this network, we see that to 
achieve this normalization we must divide the projection p: ji ® j2 ja by the square 
root of d{ji, j2, js)- The formula in Section 9.10 of the book by Kauffman and Lins 
makes it obvious that 6'(ji,j2,j3) is positive (resp. negative) when ji + j2 + js is 
even (resp. odd). This means that when we use normalized trivalent vertices, a closed 
trivalent SU(2) spin network evaluates to a real number times the product over all 
vertices v of P^+^^+ja ^ where ji, J2, J3 are the spins labelling the three edges incident to 
V. Since the spin label for each edge shows up twice in this product, it follows that the 
spin network evaluates to a real number times i"^"^ , where J is sum of the spin labels of 
all edges. 

It therefore suffices to show that whenever the 15 j symbol is nonzero, this quantity J 
is an integer. To do this, choose any vertex in the 15j symbol spin network and consider 
its images . . . ,W5 under the action of Z5 as rotations of the picture. Each edge is 
incident to exactly one of the vertices Vi, so J is the sum over 1 < i < 5 of the sum of 
spins labelling the edges incident to Vi. However, a trivalent vertex can only be nonzero 
if the sum of spins incident to it is an integer. It follows that for the 15j symbol to be 
nonzero, J must be an integer. □ 

The fact we use here about the 15 j symbols has a nice generalization: any trivalent 
SU(2) spin network evaluates to real number if its underlying graph has no loops con- 
taining an odd number of edges. The reason is that in a nonzero SU(2) spin network, 
the edges with non-integral labels can always be partitioned into disjoint loops. If these 
loops all have an even number of edges, the total of all of the spins must be an integer, 
so the above argument shows the spin network evaluates to a real number. 

Using Lemma we can easily determine the sign of the usual lOj symbol: 




Jl 
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Lemma 4. 




(-1) 



2(ji + ---+Jio) 




It follows that when 




is nonzero, it is positive (resp. negative) when ji + ■ ■ ■ + jw is an integer (resp. half- 
integer). 

Proof — If we go around this graph counterclockwise; 




each of the ten edges is incoming to exactly one of the modified Barrett-Crane inter- 
twiners. By Lemma it follows that 




(-1) 




The rest follows from Lemma |[ 



□ 
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Finally, we show that these signs cancel when we evaluate a closed spin foam of the 
sort that appears in the Barrett-Crane model. Here we assume familiarity with the basic 
definitions concerning piecewise- linear CW complexes and their homology |l^, p^ . 

Lemma 5. Let K be any 2- dimensional piecewise-linear CW complex, and let c be a 
2-cycle on K with'L2 coefficients. Let c{f) denote the value of c on any face (i.e., 2-cell) 
f , and let n[f) denote the number of vertices of the face f (which equals the number of 
its edges). Then '^fn{f)c{f) = 0. 

Proof — Write / > e when the face / has the e as one of its edges, and let (9c) (e) 
be the boundary chain of c evaluated at the edge e. Then 

E = E E ^(/) = E(^c)(^) = E = 

/ e />e e e 

□ 

We now reach our main result. Here, and in the corollary, we are using the original 
lOj symbols, so this result is independent of any choices of splittings of vertices. 

Theorem 1. Let F be a closed spin foam whose underlying 2-dimensional piecewise- 
linear CW complex K is the dual 2-skeleton of simplicial complex formed by taking a 
finite set of 4-simplices and gluing them pairwise along all their tetrahedral faces. Then 
Y\v ■^'v — 0' where v runs over the vertices of K and is the lOj symbol associated to 

V. 

Proof — We assume that the product Y[y -^v is nonzero, and show that the product 
of the signs from Lemma ^ is 1. 

Since the product is nonzero, each lOj symbol is nonzero, and therefore each Barrett- 
Crane intertwiner is nonzero. For each edge in the spin foam (corresponding to a 3- 
simplex in the simplicial complex) there are exactly four faces having it in their boundary 
(corresponding to the four triangles in the 3-simplex). Because the associated Barrett- 
Crane intertwiner is nonzero. Lemma ^ shows that the sum of the spins on these four 
faces is an integer. From the spin foam F we construct a 2-chain c with values in the 
multiplicative group {±1} by setting c(/) = (—1)^^^, where j/ is the spin labelling 
the face /. The condition above says that c is a cycle. Translating Lemma || into 
multiplicative notation, it follows that c(/)"''^^ — 1. By Lemma ^, the sign of A„ is 
X\f>vi~^)'^^^ — n/>t;'^(/)' where / runs over all the 2-cells having w as a vertex. The 
sign of n„ A, is therefore 0. <f) - D/ <^UT^^'^ = 1- □ 

In addition to a product of lOj symbols, one for every vertex, the Barrett-Crane 
amplitude for a spin foam also involves a product of factors coming from spin foam 
edges and faces. The precise formulas for these factors are somewhat controversial: 
Perez and Rovelli p2[ have proposed a modification of the original formulas given by 
Barrett and Crane |10{ . However, in both versions of the Barrett-Crane model, these 
factors are nonnegative. It follows from the above theorem that in both the original 
Barrett-Crane model and the Perez-Rovelli version, the amplitude for any closed spin 
foam is nonnegative. 

We can generalize this theorem to open spin foams as follows: 

Corollary 1. Let _F: ■0 — > ■0' be an open spin foam whose underlying 2-dimensional 
piecewise-linear CW complex is the dual 2-skeleton of a simplicial complex formed by 
taking a finite set of 4-simplices and gluing them pairwise along all their tetrahedral 
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faces except those in which the spin networks ip, ip' lie. Let G : ^/^ — > ■0' he another open 
spin foam of this type, possibly with a different underlying CW complex. If the products 
of lOj symbols over the vertices in F and G are both nonzero, then these products have 
the same sign. 

Proof — Glue F and G together along and ili' to get a closed spin foam H . Since 
the product of lOj symbols for F and G are both nonzero, the same is true of H . In 
fact, the result for H is the product of the results for F and G. By Theorem]^, the sign 
of the evaluation of H is positive, so the signs of the evaluations of F and G must be 
the same. □ 

This corollary can be sharpened using an independent argument which we now sketch. 
This sharpened version says that if the evaluation of F is nonzero, then its sign is (— 1)^'^, 
where J is the sum of the spins labelling the edges of ip ^^nd "ip' . To prove this fact, 
note that any spin foam is a composite of spin foams with only one vertex, since we can 
slide a cross-section generically from ip to ip' , so that it passes one vertex at a time. It 
thus suffices to consider the case where F contains only one vertex. For this, one can 
examine each of the possible ways in which ^ can differ from ip' when there is only a 
single vertex in F, and note that in each case, the sign of the lOj symbol at that vertex 
is (—1)^"'. For example, a portion of ip that looks like: 



This is the so-called '2-3 Pachncr move'. The sum of the ten spins involved in the lOj 
symbol is ji +^2+^3 + ^1 + ^2 + ^3 + ^ + ^711 + "I2 + ■ Since the evaluation is nonzero, 
(_l)20i+i2+j3)^ (_1)2Z (_i)2(fci+fc2+fc3) are equal. Thus the sign of the lOj symbol 
is (_i)2(i+mi-i-m2-i-m3)^ 'pjjjg jg equal to (—1)^'^, where J is the sum of the spins labelling 
ip and tp\ since most of the spins appear twice. 

There are other transitions that can occur when the cross section passes a vertex, 
such as the 0-5 and 1-4 Pachner moves. Also, the cross section can pass a cup or a cap 
in a spin foam edge, which replaces two vertices with one, or vice versa. In all of these 
cases, one can check that the signs work out as claimed. 



The Barrett-Crane model can be defined for any dimension of spacetime and any 
signature [|l6|. However, our proof of nonnegativity for spin foam amplitudes in the 




can change to: 




3. LORENTZIAN MODEL 
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4-dimensional Riemannian Barrett-Crane model relies on two miracles specific to this 
case: Spin(4) is a product of two copies of SU(2), and its simple representations factor 
as j ® j- This is what gives the weighted sum of squares in equation (p. These 
miracles do not occur in the 4-dimensional Lorentzian model, so we are unable to prove 
nonnegativity in this case. However, we have obtained some numerical evidence for it, 
which we present here. 

As in the Riemannian case, there is some controversy concerning the correct edge and 
face amplitudes in the Lorentzian Barrett-Crane model. However, in the version for 
which Crane, Perez and Rovelli have proved convergence of the partition function, 
these amplitudes are nonnegative. This is also true in all the most obvious variants. 
Given this, to prove that the amplitude of a spin foam is nonnegative, it suffices to show 
the vertex amplitudes are nonnegative. 

The vertex amplitudes in the Lorentzian Barrett-Crane model are given by the 
Lorentzian lOj symbols |ll| ]. These are obtained by evaluating Spin(3, 1) spin networks 
of this form: 




J4 



Here the edges are labelled by arbitrary nonnegative real numbers jj, which describe 
simple representations of the group Spin(3, 1) — SL(2,C). The vertices are all labelled 
with a special intertwiner called the Lorentzian Barrett-Crane intertwiner. The simple 
representation j is the space of all solutions to the equation V^/ = —{p + 1)/ on 3- 
dimensional hyperbolic space that have square-integrable boundary data on the sphere 
at infinity. Since these representations are infinite-dimensional, a careful procedure must 
be used to evaluate the lOj spin network. The end result is the following integral: 

Jl.2 / \ ^xJO.l 




J2,3 



J4.0 



(2 



1 f ^ 

-^KI / TT Kj^ ^{xy,Xw) dxidx2dx3dx4 



Here H is 3-dimensional hyperbolic space, dx stands for the volume form associated to 
the Riemannian metric on H , and Kj is the integral kernel for the projection onto the 
simple representation j, namely 

Kj{x,y) = . . 

] smh r 

where r is the distance betweeen x,y (z H . The answer is independent of the free vari- 
able xq, which may be held fixed at a convenient point in H. While this integral has 
been shown to converge 0, no 'closed form' is known, and it is difficult to compute nu- 
merically because it involves 12 variables and an oscillatory integrand. These difffculties 
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become particularly severe as the ji becorae large, since then the integrand oscillates 
more rapidly. 

We have used two methods to compute this integral numerically for certain choices of 
ji. In both methods we start by writing the integral in hyperbolic spherical coordinates 
(r, 9, (j)) G [0, oo) X [0, 2tt) x [0, tt). Then we do a coordinate transformation mapping the 
radial coordinates of the points xi, X2, x^,Xi to the interval [0, 1) via ^ ri/ {1 + ri). 
Including the Jacobian of this transformation, the result is an integral with respect to 
Lebesgue measure on a product of four copies of [0, 1) x [0,27r) x [0,7r). Finally, using 
invariance properties of the integrand, we can reduce this to a 9-dimensional integral by 
eliminating the integrals over 9i, 4>i and 62- 

The first method proceeds by a straightforward Monte Carlo calculation, uniformly 
sampling the space over which the integral is being taken, and averaging the results. 
Table 1 lists some results obtained by this method. 



jl ■■■ jio 


lOj symbol 


samples 


0000000000 


49.34 ± 1.595 


2.5 • 10'-* 


1111111111 


9.420 • 10"^ ± 3.964 ■ 10"^ 


1.0 ■ 10'-* 


2222222222 


9.138 • 10"'' ± 1.268 ■ 10"^ 


1.0 • 10** 


9999999999 


3.937 • 10"^^ ± 4.167 ■ 10"" 


1.6 • lO'-* 


19 1 1 1 1 1 1 1 1 1 


2.017 • 10"^ ± 1.498 ■ 10"^ 


9.0 ■ 10'-* 


1234512345 


1.368 ■ 10"* ± 1.814 ■ 10"" 


7.2 ■ lO'' 


1.1 1.3 0.4 2.1 1.6 1.5 0.9 1.2 1.9 0.8 


4.884 ■ 10"" ± 1.046 ■ 10"^ 


8.0 • 10^ 


1.4 15 0.1 3 2 1 1 1 0.1 8 


1.100 • 10"'' ± 1.670 ■ 10"^ 


5.2 ■ 10* 



Table 1: Lorentzian 10 j symbols — Monte Carlo results 

The second method uses the 'Vegas' algorithm. This is an adaptive algorithm that 
combines importance sampling and stratified sampling. The rather complicated details 
are explained in Section 7.8 of the book Numerical Recipes Table 2 lists some 

results obtained by this method. Each entry in the table represents a weighted average 
of independent runs of the Vegas algorithm. 



jl ■ ■ ■ Jio 


lOj symbol 


samples 


0000000000 


53.42 ± .0689 


5.0 • 10* 


1111111111 


9.526 • 10"^ ± 1.367 ■ lO"'' 


7.5 • 10* 


2222222222 


9.095 • 10"® ± 8.153 ■ lO"'' 


5.0 • 10* 


9999999999 


3.626 ■ 10"" ± 1.226 • 10"^^ 


1.3 • lO'' 


19 1 1 1 1 1 1 1 1 1 


5.602 • 10"* ± 6.934 • 10"* 


6.5 • lO'-* 


1234512345 


1.162 ■ 10"* ± 6.057 • 10"" 


1.7 ■ 10^° 


1.1 1.3 0.4 2.1 1.6 1.5 0.9 1.2 1.9 0.8 


5.029 • 10"" ± 9.258 • lO"'^ 


5.0 • 10^ 


1.4 15 0.1 3 2 1 1 1 0.1 8 


3.068 ■ 10"*' ± 6.860 ■ lO"'' 


5.5 ■ 10* 



Table 2: Lorentzian IQj symbols — Vegas results 
In both these tables, the error estimates are calculated as 
(2) A/ = 



iP) if? 



N 
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where (/) denotes the mean of the samples of the function /, N is the number of samples, 
and V is the volume of the region of integration. In the case of the Vegas algorithm 
the samples are not uniformly chosen, so the mean is a weighted one, as explained in 
Numerical Recipes. 

The reader can easily note a few things from these tables. First, the computed IQj 
symbols are always nonnegative. In two cases, namely 

(ji,...,jio) = (19,1,1,1,1,1,1,1,1,1) 

and 

(ji,...,jio) = (1.4,15,0.1,3,2,1,1,1,0.1,8) 

the computed lOj symbols are smaller than their error bars, consistent with them being 
zero or even negative. However, this is understandable, since there is a good reason 
to expect the lOj symbols to be very small in these cases. The reason is that the 
Barrett-Crane intertwiner 




becomes exponentially damped when the spins ji, . . . , leave the region defined by the 
'tetrahedron inequalities': 

3i < h+33+j4, 
h < jl+ 33+34, 

33 < 31+32+34, 

34 < 31 + 32 +33- 

This is most easily seen from the explicit formula in for the value of the Spin(4) 
spin network 



ii 




which can be thought of as the norm squared of the Barrett-Crane intertwiner. 

Second, the Vegas results generally have smaller error bars than the corresponding 
Monte Carlo results, and usually lie within the error bars of the Monte Carlo results. 
This suggests (but does not rigorously prove) that the Vegas algorithm gives more 
accurate values for the lOj symbols. 

Finally, the large sample sizes required to obtain even tolerable accuracy show the 
difficulty of these numerical computations. It will be important to develop new theo- 
retical methods to prove or disprove the conjectured nonnegativity of these Lorentzian 
IQj symbols, and also new algorithms to more efficiently compute them. 

4. Conclusions 

Given the results here together with the efficient algorithm for computing the Rie- 
mannian lOj symbols described in a companion paper [ p^ , it is now practical to study 
the Riemannian Barrett-Crane model using the Metropolis algorithm. Some results 
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along these lines will appear in another paper For the Lorentzian model we still 
need better formulas for the lOj symbols, both for proving nonnegativity and for nu- 
merical computations. 

The physical significance of Corollary 1 is still somewhat mysterious. Multiplying the 
usual spin network states by suitable signs, it means we can find orthonormal bases of 
the kinematical Hilbert spaces of any triangulated compact 3-manifolds S and 5" such 
that for any triangulated cobordism Af : S S', the matrix elements of the operator 
Z{M) : Z{S) Z{S') are always a sum of nonnegative terms, one for each spin foam 
F: ip ^ ij}' iu the dual 2-skeleton of M. At first it was suggested that this means there is 
no destructive interference between spin foams in the Riemannian Barrett-Crane model. 
Subsequent discussions have clarified this issue somewhat. Indeed, if a unitary operator 
has all nonnegative matrix elements in some basis it must be rather trivial, simply acting 
as a permutation of the basis. However, thanks to the 'problem of time' in quantum 
gravity, it is more natural to use a sum over spin foams to compute, not a time evolution 
operator, but the projection from the space spanned by spin networks onto the physical 
Hilbert space of quantum gravity |2^. There is no reason for this to be unitary; it 
should be something more like a positive operator. There is thus no contradiction if 
it has nonnegative matrix elements. Indeed, this also occurs in other theories, such as 
the quantum mechanics of a relativistic free particle, and also quantum gravity in 2 -I- 1 
dimensions. However, we still lack a general principle explaining why all spin foams 
F: ip ^ should have real amplitudes of the same sign. 
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